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Mechanical stretching of secondary structures is studied 
through molecular dynamics simulations of a Go-like model. 
Force vs. displacement curves are studied as a function of the 
stiffness and velocity of the pulling device. The succession of 
stretching events, as measured by the order in which contacts 
are ruptured, is compared to the sequencing of events during 
thermal folding and unfolding. Opposite cross- correlations 
are found for an a-helix and a /3-hairpin structure. In a 
tandem of two a-helices, the two constituent helices unravel 
nearly simultaneously. A simple condition for simultaneous 
vs. sequential unraveling of repeat units is presented. 



INTRODUCTION 

Weak noncovalent bonding forces govern functioning 
and structural cohesion in cells. Direct measurements of 
these forces through mechanical means has recently be- 
come an important tool in studies of biological molecules. 
There is a variety of techniques for probing forces in the 
pico- and nano-Newton range (lj such as atomic force 
microscopyj |-|| , optical tweezers ||J^], the surface force 
apparatus |B|, micropipette aspiration and the quartz 
microbalance ||lfj| , p4 |. As examples of recent achieve- 
ments, we may list elucidation of the nature of interac- 
tions of a chaperone protein (HIV-1) with DNA through 
stretching of a strand of the DNA with optical tweezers 
|l2f and discovery of stick-slip motion when two strands 
of a DNA double helix are pulled apart Jl3| ]. 

The techniques used in mechanical unfolding of in- 
dividual biological molecules rely on tethering of the 
molecule between movable surfaces. This tethering is 
relatively easy to accomplish with long molecules such as 
DNA and giant proteins such as titin which are nat- 
urally built as a tandem array of many globular domains. 
For shorter molecules, the pulling surfaces interact and 
affect the pulled molecule in a way that makes the data 
hard to interpret. In order to extend the method to single 
domained proteins, Yang et al. E3] have recently devel- 
oped a method of synthesizing identical repeats of protein 
molecules in the solid state which were then studied us- 
ing a modified scanning force microscope. This technique 
has been applied to T4 lysozyme. 



At this moment, experimental data on the mechanical 
unfolding of the secondary structures of proteins are not 
available. However, data on periodically repeated pro- 
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teins and even individual proteins may become available 
in the near future. From a theoretical point of view, it 
is important to gain an understanding of the basic un- 
folding mechanisms of simple structures and to develop 
analytical tools that could then be used for large pro- 
teins. This process is facilitated by considering simple 
models that allow a rapid exploration of parameter space. 
Our choice in this paper is to analyze Go-like models 
which emphasize the importance of native conformations 
and treat non-native interactions only schematically. The 
Go-like models |l5j , though coarse-grained, are fairly re- 
alistic 16 1 in their kinetic properties and allow for a 



thorough characterization and comparison of mechani- 
cal, equilibrium and folding properties in a straightfor- 
ward manner. This kind of full characterization is diffi- 
cult to achieve in all-atom models with the Amber Jl7|] 
or CHARMM force fields which nevertheless are well 
suited to studies of mechanical stretching. 



cific contacts are broken. Both the force-displacement 
curve and the order of unfolding events depend on the 
stiffness and velocity of the cantilever. We next discuss 
studies of folding, where temperature plays an essential 
role. The sequencing of folding events depends on T, and 
smooth and simple pathways are only found near an opti- 
mal temperature denoted by T min p^-|2?[| . The sequenc- 
ing of folding events near T m i„ is contrasted with that 
of stretching events for different protein structures. We 
find that both sequencings are governed primarily by the 
contact order ]2g|-|30[] , i.e. by the distance between two 
amino acids along the sequence of the protein. However, 
the cross-correlations between thermal and mechanical 
sequencings are opposite for the two simple cases consid- 
ered: H16 and B16. Only in the latter case do folding 
and stretching occur in the opposite order, as envisioned 
by, e.g., Bryant et al. GjJ. In general, the thermal and 
mechanical pathways can be very different. 



The idea that mechanical unfolding experiments on 
proteins have the potential to provide insights into the 
relevant folding pathways is what motivated Bryant et 
al. jl9) to carry out all-atom (CHARMM-based) simu- 
lations of the C-terminal hairpin of protein G, the fold- 
ing of which has been previously studied experimentally 
by Munoz et al. pO| , pl| . They have found that, under 
low pulling forces, breakdown of hydrogen bonds pre- 
cedes dissociation of the hydrophobic cluster. Their in- 
terpretation of this finding is that thermal folding should 
proceed in the opposite order to mechanical unfolding. If 
so, then the zippering folding mechanism | pl| would be 
less favored than one in which a hydrophobic cluster is 
formed first. This prediction remains to be tested. 



Another quantity that we study here is what we pro- 
pose to call an irreversibility length, Li r . If one stud- 
ies folding from a fully extended conformation, then 
one finds that the characteristic folding time diverges as 
T — >0. Thus a fully stretched protein will not fold back 
to the native state at low temperatures. On the other 
hand, a protein that is pulled only slightly will return 
to its native shape on release. There must then be a 
characteristic stretched length of the protein, Li r , which 
separates the two behaviors. We demonstrate that Li r 
does indeed exist and find that it is substantially larger 
for B16 than for H16. Furthermore, the folding time for 
lengths less than Li r is a complicated function of the me- 
chanical extension. 



Here, we explore properties of Go models of proteins 
through molecular dynamics simulations. We consider 
the variant in which contact interactions are described 
by Lennard-Jones potentials. The simulations include a 
Langevin noise term which both mimics presence of a sol- 
vent and controls the temperature, T. This paper focuses 
on an a - helix of 16 monomers, denoted as H16; a (3 - 
hairpin of 16 monomers, B16; and a double repeat of the 
a - helix, H16-2. The companion paper |^2| describes a 
similar analysis for titin. 

We first study mechanical unfolding at nearly zero tem- 
perature. This choice of T minimizes fluctuations and 
rate dependence, and most simply reveals the effects of 
the structure of the energy landscape. The results should 
be equivalent || to fast stretching at higher T, and tem- 
perature dependence will be considered in subsequent 
work. The protein is stretched by a Hookean cantilever 
and the force is plotted as a function of the cantilever dis- 
placement. We also characterize the stretching process 
by studying the succession of unfolding events, which are 
described by the cantilever displacements at which spe- 



We also consider a tandem arrangement of two a- 
helices and find that the constituent helices unravel al- 
most simultaneously whereas in titin | ^2f the unraveling 
is serial in nature. Simple criteria for the two types of 
behavior are described. 



MODEL AND METHOD 

The model we use is described in detail in references 
p5|- p7| , j3T|| . For simplicity, we consider the variant where 
steric constraints associated with dihedral and other an- 
gles are ignored. Briefly, a protein is modeled by a chain 
of identical beads which correspond to the locations of 
the Ca atoms. The consecutive beads interact through 
the potential |32| 
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where r^j+i = |r^ — is the distance between two con- 
secutive beads, d — 3.8 A, fci = e/A 2 k 2 — 100e/A 2 , and 
e is the characteristic energy parameter corresponding to 
a native contact. The anharmonic term in equation (1) 
prevents energy localization in specific phonons and thus 
accelerates equilibration p^ j. 

The interaction that governs the native contacts (de- 
fined as those Ca that are not immediate neighbors, but 
are no further than 7. 5 A apart in the native structure) is 
chosen to be of the Lennard- Jones type (see e.g. pjJ) 
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The parameters <xy are chosen so that each contact in 
the native structure is stabilized at the minimum of the 
potential, and a = 5A is a typical value. As a technical 
criterion for determining when a native contact forms or 
breaks during the time evolution, we adopted the cutoff 
value of 1.5c7jj. The non-native contacts are described by 
purely repulsive potentials. These are obtained by eval- 
uating V NAT with a length parameter a, truncating the 
potential at its minimum (2 1 / 6 cr), and shifting it to have 
zero value at this cutoff distance. 

Figure 1 illustrates the forms of the potentials for 
the a-helix. When studying the folding times, we have 
adopted a simplified approach in which a protein is con- 
sidered folded if all beads that form a native contact are 
within the cutoff distance of 1.5(7^ instead of making a 
more precise delineation of the native basin as in ref . p6[ . 
This will allow for a more meaningful comparison with 
the results on titin p2[ . 

The beads are coupled to Langevin noise and damp- 
ing terms to mimic the effect of the surrounding solution 
and maintain constant temperature T. The equations of 
motion for each bead are 



- 7 r + F c + T 



(3) 



where m is the mass of the amino acids represented by 
each bead, F c is the net force due to the molecular po- 
tentials and external forces, 7 is the damping constant, 
and r is a Gaussian noise term with dispersion v / 2 7 /cb 7 T. 
We measure time in units of the characteristic period 
of undamped oscillations in the Lennard-Jones potential 
t = yj mcr 2 /e. Using typical values for the average amino 
acid mass, length and binding energy yields 3ps as an 
estimate of r. According to Veitshans ct al. fl34| , realis- 
tic estimates of damping by the solution correspond to 
a value of 7 near 50 m/r. However, the folding times 
have been found to depend on 7 in a simple linear fash- 
ion for 7 > m/r p5| , p6| , p5[ . Thus in order to accelerate 
the simulations, we work with 7 = 2m/r. The equa- 
tions of motion are solved by means of the fifth order 



Gear predictor-corrector algorithm |36| with a time step 
of 0.005r. 

In order to pull the protein apart, we attach both of 
its ends to purely harmonic springs of spring constant k. 
We focus on three cases: a) the stiff spring: k > 60e/A 2 
b) the soft spring: k = 0.12e/A 2 , and c)the very soft 
spring: k = 0.04e/A 2 . The outer end of one spring is 
held stationary, and the other is pulled at a fixed rate 
v p . This models stretching by a Hookean cantilever with 
stiffness fc/2, since the two springs add in series. Wc 
also performed simulations at constant force, which cor- 
responds to the limit of infinitely weak springs. However, 
the unwinding of the proteins occurs in an " all or noth- 
ing" fashion in this limit, and little information can be 
extracted. 

There are many ways to pick the pulling direction, and 
all of our mechanical results correspond to a case in which 
the extension is implemented parallel to the initial end- 
to-end vector. In most cases, we pull the spring very 
slowly - at a constant rate of v p — 0.005A/r. There is 
actually very little dependence of the results on pulling 
rate until one considers rapid rates. For instance, in- 
creasing v p by a factor of 50 produces almost no change 
in the force. Substantial rate dependence begins when 
v p is increased by a factor of 100 to 0.5 A/r, and this 
case is denoted as a "fast" stretch in the following sec- 
tion. The instantaneous pulling force F is the extension 
of the pulling spring times the spring constant k. Plotted 
values of F are averaged over It. The standard pulling 
velocity is low enough that the force equilibrates along 
the chain and almost the same force is obtained from the 
extension of the spring whose end is fixed. Drag terms 
lead to a significant difference in these forces at higher 
velocities. The force is plotted versus the cantilever dis- 
placement d — Vpt, where t is the total pulling time. 



RESULTS AND DISCUSSION 

a-helix 

Figure 2 illustrates the process of mechanical unfolding 
for H16. It clearly shows that unfolding starts at both 
ends and then proceeds to the center. This is precisely 
the ordering of events during thermal folding [p5| and 
not the inverse of this ordering as seen for the C-terminal 
hairpin of protein G [ pi| . However, the underlying rea- 
sons for the observed ordering during folding and unfold- 
ing are different. Folding starts at the ends because they 
diffuse more rapidly and are thus more likely to fall into 
a contact situation, while unfolding starts at the ends 
because there are fewer binding forces there. 
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The force vs. cantilever displacement curves are shown 
in Figure 3 for slow (solid lines) and fast (dashed lines) 
displacement rates. The curves are truncated when the 
helix is fully stretched, and any further displacement re- 
sults in rapid growth in F followed by rupture of the pro- 
tein backbone. The upper (lower) panels are for the stiff 
(soft) pulling springs. In both cases unfolding produces 
a sequence of stick-slip events. The force rises linearly 
while the protein is trapped in a given local energy min- 
imum, and then drops rapidly as one or more contacts 
breaks. The slope of the upward rise is the combined 
stiffness k to t of the protein k p and cantilever k/2. Since 
the two are in series, k^ ot — k^ 1 + 2fc _1 . In the soft 
spring case the cantilever dominates and the slope of the 
upward ramps is k/2. For the stiff cantilever case the 
internal stiffness of the protein dominates. Variations 
within and between local minima lead to changes in the 
slope of the ramps, with k p varying between about 0.3 
and 0.6 e/A 2 . Once k is larger than these values it has 
little influence on the curves. 

Each upward ramp ends when one or more contacts 
break. The force drops sharply until the protein reaches 
a new metastable state and a new upward ramp begins. 
In the low velocity case (solid lines), v p is much lower 
than the velocities produced by contact breaking, and 
rupture occurs at a nearly constant cantilever position. 
In this limit, the force drop is roughly equal to ktot times 
the change in protein length during the jump between 
metastable states. For a stiff cantilever (top panel), the 
failure of each contact produces a large drop in the force. 
The first two peaks correspond to breaking of the two end 
contacts. The force is lower than for later events because 
the ends have fewer native contacts. Rupturing of the 
next series of bonds proceeds in an essentially periodic 
pattern because each ruptured bond has the same envi- 
ronment. When the remaining helical segment is short 
enough, failure affects bonds across its entire length, lead- 
ing to two higher peaks. In the final stages (i.e. d > 28A 
), all the coils have been broken and the series of small 
force peaks is due to breaking of higher order bonds. 

When a soft spring is used, the drop in force due to 
each event is smaller. If the threshold force for an event 
is lower than that for the previous event, the force may 
not drop below this threshold. This can cause several 
bond ruptures to accumulate into a single orchestrated 
event. The low velocity curve in the lower panel of Fig. 
3 has the same initial sequence of peaks as the upper 
panel: Two small peaks are followed by several at the 
same higher force. However, those later peaks that are 
well below preceding peaks in the top panel are absent 
in the bottom panel. The strength of the contacts bro- 
ken in these stages would be difficult to extract if a soft 
cantilever were used. 

When the pulling velocity is comparable to the rapid 
motions produced by bond rupture, the cantilever mo- 



tion can produce a substantial change in force during an 
unfolding event. This can also cause events to accumu- 
late as shown in both panels. The increase in speed also 
produces a larger drag force from the surrounding solu- 
tion (represented by the Langevin damping). This shifts 
the force curves to higher values. 

The optimal temperature for folding of H16 has been 
established to be T min = 0.3e/fcs ^5|. The sequenc- 
ing of thermal folding and unfolding events at T m i n is 
shown in Figures 4 and 5. The former Figure considers 
establishment of the contacts of the i, i + 4 type, i.e. the 
hydrogen-bonded contacts, whereas the latter is for the 
i, i + 3 contacts. The time for establishing of a given 
contact is denoted by t c . These times are symmetrically 
arranged around the center of the helix and are shortest 
at the ends. We have also determined times for thermal 
unfolding, t u , defined as times at which the contact is 
gone for the first time. Values of t u in Figures 4 and 
5 are averaged over 1500 different trajectories which all 
start in the native state. We notice that t u is an or- 
der of magnitude longer than t c but it is also arranged 
symmetrically around the center of the helix. 

Figures 4 and 5 also show the displacements, d u , at 
which each bond ruptures during mechanical unfolding 
for stiff (closed circles) and soft (open circles) cantilevers. 
These curves do not have the same symmetry as the ther- 
mal folding and unfolding curves. As noted above, the 
end bonds break first because they have fewer native con- 
tacts. Subsequent bonds have the same number of con- 
tacts and should break at the same force. However, the 
bonds near the pulling end (large i) tend to break first 
due to the presence of a small extra drag force. This is 
independent of the nature of the cantilever except that 
the soft spring yields uniformly larger d u at which a bond 
breaks. 

Despite the lack of symmetry in the mechanical data, 
the folding times and contact breaking distances are 
clearly correlated. This is shown in Figures 6 and 7 for 
the stiff and soft spring respectively. In each figure, con- 
tacts breaking at later times tend to break at larger dis- 
placements. 

Two helices in tandem 

We now consider two H16 helices connected in series 
by one extra peptide bond and stretched from one end. 
Figure 8 shows a snapshot of a partially unfolded tandem 
conformation. It indicates that the two helices unfold si- 
multaneously with some phase shift between them. This 
is also seen in the F — d curves shown in Figure 9 where 
the stick-slip patterns essentially double each feature seen 
in Figure 3. This behavior is quite distinct from what 
happens when stretching titin where the domains unfold 
one at a time p3| . The basic reason that the helices un- 
fold simultaneously is that the force to break contacts 
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rises smoothly during the unfolding process. The heights 
of the force peaks only drop in the very late stages of 
growth when the coils are all gone. In the case of titin, 
one of the early peaks is higher than subsequent peaks. 
Once this contact breaks in one of the repeat units, there 
is a series of weaker bonds that can continue to rupture 
within that unit. These contact failures keep the force 
from rising back to a level that would initiate failure of 
the strong bonds in other repeat units. 

The simultaneous unwinding of the two helices is also 
seen in Figure 10 which is an analog of Figure 4 for the 
single helix (minus the data on thermal unfolding). The 
distance for contact rupture (for contacts) through 

stretching shows two skewed peaks, each centered in the 
vicinity of the centers of the individual helices. In con- 
trast, the average times for unfolding at T — 0.3e/fce are 
peaked not at the centers of the helices but at the very 
center of the whole system, i.e. around the peptide bond 
that connects the helices. Thus the simple correlation 
between d u and t c that was seen in Figures 6 and 7 is 
lost. Instead one finds a two-legged correlation that is 
shown in Figure 11. Note also that all of the contacts 
(all are short ranged and are grouped into three types: 
i, i + 2, i, z+3, and i, z+4) break throughout the full range 
of the displacement of the cantilever. Some bonds of a 
given kind break early, some break late. We shall see in 
the companion paper [ p2[ that failure of long range bonds 
shows a definite correlation with the displacement. 

/3-hairpin 

The stretching of the /3-hairpin B16, shown in Figure 
12, consists of a gradual removal of the "rungs" of the 
"ladder" that form the hairpin, starting from the free 
ends. Physically, these rungs represent hydrogen bonds 
and they correspond to contacts 1-16, 2-15, 3-14, 7- 
10. There are other contact forces in our Go model and 
they provide further stabilization of the structure. These 
other bonds bind bead 1 with bead 15, bead 2 with 14 
and 16, etc. 

Plots of F versus d during unfolding at low v p are 
shown in Figure 13. All show regular stick-slip fea- 
tures. In this respect, our results are very similar to 
those obtained by Bryant et al. |l^] with full atom sim- 
ulations. Thus our simplified model reproduces the fea- 
tures present in the more realistic calculation. Further- 
more, since our model incorporates the native conforma- 
tion but not the hydrophobic or polar properties of the 
amino acids, we suggest that the latter are not explic- 
itly crucial in the mechanical unfolding of the hairpin. 
The stiff and soft springs produce the same sequence of 
stick-slip peaks, but the slope of the ramps and depth of 
the drops are smaller for the soft spring. After the first 
peak, peaks come in pairs where the second peak has a 
lower height. When a very soft spring is used, these pairs 



merge into single large events as described above. 

The folding properties of B16 are illustrated in Figure 
14. This system has been studied in detail in ref. p5| , 
where the native basin has been accurately determined 
through a "shape distortion technique" |37j which pro- 
duces T min of order 0.07e//c_B- If the folding criterion 
is based on just establishing the native contacts, then, 
in the case of B16, there is a very broad dependence of 
the folding time on temperature and the kinetics of fold- 
ing at 0.07e/ks is almost the same as at, say, 0.3 e/fc^. 
Nevertheless we study the system at the previously de- 
termined T m i n . Note that even with the contact-based 
criterion for folding, the folding time for B16 is still con- 
siderably longer than for H16. 

Figure 14 shows that the sequencing of folding events 
in B16 is exactly opposite to the succession of contact 
breakage upon stretching: B16 starts folding from the 
turn (the result that has been found both experimentally 
f20f and theoretically [^,|8|) whereas both mechanical 
and thermal unfolding start at the free ends. Thus, in 
contrast to the a-helix, the mechanical unfolding of the 
/3-hairpin is the inverse of the folding process. 

Figure 15 shows d u as a function of the time needed 
to establish the contact during folding. Here, in addition 
to the "rung" contacts, the remaining contacts are also 
shown. Since contacts rupture at a fixed force, the soft 
spring data are shifted to larger displacements than the 
the stiff spring data. However both sets of data show 
a clear anticorrelation between thermal folding and me- 
chanical unfolding that is in a sharp contrast to the re- 
sults for the a-helix. 

Irreversibility length 

We now consider pulling of a protein at a constant 
slow rate and then releasing it. We ask what is the 
time required to fold back to the native state at T = 0. 
There must be a limit to the extension beyond which the 
protein misfolds on release. Figure 16 shows that this 
limit indeed exists. The dependence on cantilever stiff- 
ness is minimized by plotting the refolding times against 
the end-to-end distance L of the protein rather than the 
cantilever displacement. For both stiffnesses the refold- 
ing times are found to be non-monotonic functions of 
L. We interpret this as being due to inertial effects. The 
more stretched the protein is with a given set of contacts, 
the more potential energy is available. When the protein 
is released, the energy is converted into kinetic energy 
that speeds the contraction of the protein and aids it in 
getting over subsequent energy barriers. 

We identify the irreversibility length L; r with the max- 
imum value of L where refolding occurs. For H16, L; r is 
about 37A, or 1.6 times the native state end-to-end dis- 
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tance of 22.62 A. The change in length is 14. 4A which is 
very close to the displacement of the stiff cantilever at the 
onset of irreversibility di r = 14. 9A . The displacement of 
the soft cantilever, di r — 37. 4A , is larger because the 
cantilever stretches more in order to apply enough force 
to reach L ir . Examining Figure 3, we see that both val- 
ues of di r correspond to the displacement after the sixth 
peak in the respective force curve. Thus the same set 
of broken bonds is required to produce irreversibility for 
either cantilever stiffness. 

For B16, the native L is only 5 A and the stretching 
factor to Li r is substantially larger, around 11.6. The val- 
ues of d^ for stiff and soft cantilevers are d ir — 52.9 and 
65. 9A, respectively. From Figure 13 we see that in both 
cases the irreversibility point is just past the last peak 
in the force curve. Since the protein is fully stretched 
at this point, any native contacts are enough to ensure 
refolding. 

The misfolded conformations, that are obtained on re- 
folding beyond the threshold, are shown in Figure 17. In 
the case of B16, the turn region freezes into the wrong 
configuration which is almost straight. In the case of 
H16, the first turn coils with the wrong chirality. 



CONCLUSIONS 

We have studied the force - displacement curves for sec- 
ondary structures of proteins for two models of cantilever 
stiffness and several pulling speeds. A series of stick-slip 
events is observed as contacts break. Stiff cantilevers 
pulled at low rates provide the most detailed informa- 
tion about the breaking of individual contacts. Multiple 
contact ruptures merge into single events when the stiff- 
ness is decreased or the speed is increased. 

The simple expectation that mechanical unraveling 
should proceed in the inverse order from thermal fold- 
ing is only confirmed in the case of the /3-hairpin. In 
the case of the a-helix, unraveling and folding follow the 
same order. When multiple helices are connected in tan- 
dem, the correlation becomes even more complex. The 
two helices unravel simultaneously with each helix uncoil- 
ing from both of its ends. In contrast, folding occurs first 
at the outer ends of the pair of helices. In general there 
is no reason to expect a simple correlation between ther- 
mal folding and mechanical unfolding of proteins. In the 
following paper we examine similar issues of mechanical- 
thermal correlations for a protein with a significant num- 
ber of long-ranged contacts. 
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FIG. 1. The potentials used to construct a Go model 
of the a-helix H16. The interactions are between the beads i 
and i + k. For fc=l this is the anharmonic tethering potential. 
The contact corresponding to fc=5 is non-native and is thus 
purely repulsive. The remaining contact interactions are of 
the Lennard- Jones form. 
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FIG. 2. Snapshot pictures of stretching of the a-helix H16. 
The left end is anchored elastically and the right end is pulled 
by the stiff cantilever. The numbers indicate the cantilever 
displacements. In the upper figure the helix is still almost 
fully folded. 
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FIG. 4. Sequencing of events as measured through forming 
or breaking bonds of the i,i + 4 kind. The x-axis shows i - the 
monomer number along the chain. The moving cantilever is 
attached to i=16. The y-axis shows t c , t u , and d u - The first 
quantity is the mean time needed to establish the contact at 
Tmin (based on 1000 trajectories). The second quantity is the 
first unfolding time at T m i n for a bond under the conditions 
of no pulling force (based on 1500 trajectories). The third 
quantity, denoted by the circles, is the displacement where 
the contact is broken during mechanical unfolding at T—0. 
The solid (open) circles correspond to slow pulling by the stiff 
(soft) spring. Values of d u for the soft cantilever are divided 
by 2. Overall, the size of the symbol is a measure of the error 
bars, and all lines are guides to the eye. 
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FIG. 3. Force vs. displacement for H16 with stiff and soft 
cantilevers. The solid (dotted) lines correspond to the slow 
(fast) pulling rates. 
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FIG. 6. Stretching distances at which a bond rupture 
takes place (from Figures 4 and 5) plotted vs. average time 
needed to establish contact on folding. This is the case of a 
stiff spring which is being pulled slowly. 



FIG. 8. Conformation of two a-helices connected in series 
after moving the cantilever by the distance indicated. 
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FIG. 9. Force vs. displacement for H16-2 - two a-helices 
H16 connected in series and pulled slowly. 
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FIG. 10. Sequencing of events in H16-2 as measured 
through the bonds of the i, i + 4 kind. The symbols have 
meanings as in Figure 4. 1000 trajectories were used in the 
studies of folding. 
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FIG. 11. Stretching distances at which a bond rupture 
in H16-2 takes place vs. average time needed to establish 
contact on folding. The cantilever is stiff and pulled slowly. 
The contacts corresponding to the i, i + 2 type are shown 
by stars. The remaining symbols differentiate between the 
contacts present in the first helix (circles) and those present 
in the second helix (squares) - the one which is closer to the 
cantilever. The open circles and squares correspond to the 
contacts of the i, i + 3 type and the filled circles and squares 
to i, i + 4. 
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FIG. 12. Snapshot pictures of stretching of the /3-hairpin 
B16 for d equal to 5 and 20 A. 
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FIG. 13. 



The force 



displacement curves for the 



/3-hairpin B16 obtained at slow pulling velocities. 
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FIG. 15. Stretching distances at which a bond rupture 
in B16 takes place vs. average time needed to establish con- 
tact on folding. The cantilever is pulled slowly. The solid 
circles correspond to the "rung" contacts. The open circles 
correspond to the remaining contacts. There are degeneracies 
related to these other contacts. For instance, 2-14 forms and 
breaks at essentially the same moment (statistically) as 3-15. 
Examples of other such pairs are 3-13 with 4-14 and 6-12 with 
5-11. 



FIG. 14. Similar to Figure 4 but for the hydrogen rung-like 
bonds in B16. The bonds are identified by the index i that 
they connect to. The thermal data are based on 1000 trajecto- 
ries and are collected at T m in = 0.07e/fcs. The flat character 
of the data corresponding to thermal unfolding is expected to 
turn into a steeper dependence at higher temperatures. 
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FIG. 16. Refolding times after stretching to the indicated 
end-to-end distance. The main figure is for B16 and the inset 
for H16. The solid (dotted) lines correspond to pulling by a 
stiff (soft) cantilever. The curves end at Li r . To the right 
of the data points shown, the protein does not return to its 
native state. For B16, the corresponding threshold values of 
the tip displacement, dj r , are equal to 52. 9A and 65. 9A for the 
stiff and soft cases respectively. For H16, the values of di r for 
the stiff and soft cantilevers are 14. 9A and 37. 4 A, respectively. 
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FIG. 17. Conformations corresponding to the misfolded 
proteins after stretching just beyond the irreversibility thresh- 
old. 
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